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This is a summary of the proof by G.E. Coxson pQ that P-matrix recognition is co-NP- 
complete. The result follows by a reduction from the MAX CUT problem using results 
of S. Poljak and J. Rohn [5]. 

1 Considered problems 

Our main interest is the complexity of deciding whether an input matrix is a P-matrix. 
A P-matrix is a square matrix M E M. nxn such that all its principal minors are positive. 
Such matrices were first studied by Fiedler and Ptak [2]. 

P-MATRIX 

Instance: A square matrix M G Q nxn . 

Question: Are all the principal minors of M positive? 

To start with, we use a well-known combinatorial problem. 

SIMPLE MAX CUT 

Instance: A graph G = (V, E), a positive integer K. 

Question: Is there a partition of the vertex set V into sets V\ and V2 such 
that the number of edges with one end in V\ and the other end 
in V% is at least K? 

Garey, Johnson and Stockmeyer [5] showed that the SIMPLE MAX CUT problem is 
NP-complete. 

The reduction from SIMPLE MAX CUT to P-MATRIX uses two intermediate steps. 
The first of them is the computation of the r-norm of a matrix. 
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For an arbitrary matrix A G M nxn , let 

r(A) = max{z T Ay : z,y G {-1, 1}"} . 
Remark. The function r is a matrix norm. 

Proof. For an arbitrary square matrix A, we have r(A) > because z T Ay = — (—z) T Ay. 
Moreover if r(A) = 0, then z Ay = for all choices of z,y G {—1, l} n , hence A = 0. If 
fceR, then z T (kA)y = k ■ z T Ay, so r(kA) = \k\ ■ r(A). 
Let A,B e M. nxn . Then 

r{A + B) = max{z T (A + B)y :y,z£ {-1, l} n } = max{z T Ay + z T By :y,z£ {-1, l} n } 
< max{z T Ay : y,z G {-1, 1}™} + max{z T By : y,z G {-1, 1}™} 

= r(A)+r(B). 

Thus r is also subadditive. □ 
The decision problem corresponding to r-norm computation is defined as follows. 

MATRIX R-NORM 

Instance: A matrix A £ <Q nxn and a rational number K. 
Question: Is r(A) > Kl 

For the last of the decision problems considered here, we need the notion of matrix 
interval. If A_ and A + are nxn real matrices such that A- < A + (that is, for each 
r and s we have (^4_) rjS < (A + ) r ^ s ), then the matrix interval [A-,A+] is the set of all 
matrices A satisfying A_ < A < A + . 

A matrix interval is singular if it contains a singular matrix; otherwise it is non- 
singular. 

The decision problem we consider consists in testing whether a given matrix interval 
is singular. We will see that this is a computationally hard problem even when the 
difference A + — A- has rank 1. 

RK1-MATRIX-INTERVAL SINGULARITY 

Instance: A non-singular matrix A £ Q nxn and a non-negative matrix A G 

Qnxn of rank L 

Question: Is the matrix interval L4 — A, A + A] singular? 

The rest of this exposition contains three polynomial reductions of these problems, 
ultimately proving that P-MATRIX is co-NP-complete. 



'This object is usually called an interval matrix. Since it is actually an interval and not a matrix, I 
beg the reader to pardon my decision to call it an uncommon but appropriate name. 
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2 Reduction from SIMPLE MAX CUT to MATRIX R-NORM 



Let G = {V,E) be an undirected graph with n = \V\ and let £ = 2\E\ + 1. If A(G) is 
the adjacency matrix of G, define A = I ■ I n — A{G). Thus 

£ if u = v, 

-l ituveE, 

otherwise. 

Observe that for y,z E {— 1, l} n we have z T Ay < y T Ay because of the choice of £. 
Hence r(A) = y T Ay for some y € {— 1, l} n . 

Let S C V be defined by S = {u : y u = 1} and let w! be the number of edges of G 
with one end in S and the other end in V \ S. In this way, m! is the size of the cut 
defined by S and V \ S. 

Then 

y T Ay = n£ + 4m' - 2|^7| 

and therefore there is a cut in G of size at least K if and only if r(A) > n£ — 1\E\ + 4K. 

The described reduction (by Poljak and Rohn [5]) establishes the hardness of comput- 
ing the r-norm. 

Theorem 1. MATRIX R-NORM is NP-complete, even if input is restricted to non- 
singular matrices. 

Proof. It follows from the reduction above that MATRIX R-NORM is NP-hard. Observe 
that by the choice of £ the matrix A in the reduction is strictly diagonally dominant and 
thus non-singular. 

A non-deterministic Turing machine can guess the values of y,z € {— l,l} n and check 
in polynomial time that z T Ay > K, so the problem is in the class NP. □ 

3 Reduction from MATRIX R-NORM to 
RK1-MATRIX-INTERVAL SINGULARITY 

For a matrix A € M nxn define 

Po(A) = max{|A| : A is a real eigenvalue of A} 

and set po(A) = if A has no real eigenvalue. 

Further for a vector y € M. n define D{y) to be the diagonal n x n matrix with diagonal 
vector y. 

The following fact was proved by Rohn [6]. 

Lemma 2. Let A be a real non-singular n x n matrix and let A be a real non-negative 
n x n matrix. Then the matrix interval [A — A, A + A] is singular if and only if 
Po(A- 1 D(y)AD(z)) > 1 for some y,z£ {-1,1}". 
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Proof. For y, z G { — 1, l} n let A^ denote the matrix D(y)AD(z). 

First suppose that A^ 1 A y>z has a real eigenvalue A such that |A| > 1 and A~ l A y ^ z x = 
Ax for some y,z G { — 1,1}™ and a non-zero vector x. Then 

(1 - ^A^Ay,,) x = 0, 

(A-lA^)x = 0. 

Hence A — (l/A)A 2/;2 is a singular matrix in the interval [A — A, A + A] because 

\\Ay, g \ = \\D(y)AD(z)\<A. 

Therefore the interval [A — A, A + A] is singular. 

To prove the converse, suppose that B is a singular matrix, B G [A — A, A + A]. Let 
x be a non-zero vector for which Bx = 0. 

For i = 1, 2, . . . ,n set 

_ 

* (A|s|V 

We claim that t G [0, l] n . Indeed, |Ac| = |(v4 - B)x\ < A\x\ because Bx = and 
Be[A-A,A + A]. 

Moreover, set z = sgnx. Then D(z)x = \x\ and 

(A - A t:Z )x = Ax- D{t)AD{z)x = Ax - D{t)A\x\ = 

by the definition of t. Thus the matrix A — At )Z is a singular matrix in the interval 
[A- A,A + A]. 

Define ip(s) = det(A— A S)Z ). The function ip is affine in each of the variables s±, . . . , s n . 
Since ip(t) = det(A—A t>z ) = 0, either there exists y G {—1, l} n such that det(A— A y ^ z ) = 
0, or there exist y,y' G {-1, 1}™ such that det(^4 - A y>z ) ■ det(A - A y ,^ z ) < 0. 

In the latter case, without loss of generality we may assume that det A-det(A — A y%z ) < 
0. The function <fi defined by 4>(a) = det(^4 — aA y ^ z ) is continuous and <^>(O)0(1) < 0, so 
<j) has a root in (0, 1). 

In either case, there exist y G {—1, l} ra and a G (0, 1] such that det (A — aA y>z ) = 0. 
Then 

det (±A - A y:Z ) = 0, 
det (i J - A- x A y , z ) = 0, 

hence ^ is a real eigenvalue of the matrix A~ l D{y)AD{z) and ^ > 1, as we were 
supposed to prove. □ 

This lemma provides a useful connection between singularity of matrix intervals and 
a parameter po dependent on the two matrices A, A that define the interval. Next we 
establish a connection between po and the r-norm of matrices. 

From now on let 1 be the all-one vector (1, 1, . . . , 1) G M. n and let J = 1 • 1 T be the 
all-one n x n matrix. 
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Lemma 3. Let A G IR nxn be a non-singular matrix, let a be a positive real number and 
let A = a J. Then 

max{po(AD(y)AD(z)) :y,ze {-1, l} n } = a • r(A). 

Proof. First observe that D(y)AD(z) = a • D(y)l ■ l T D(z) = a ■ yz for arbitrary 
y, z G {— 1, l} n . If A is a non-zero real eigenvalue of a ■ Ayz T and x is a non-zero vector 
such that 

a ■ Ayz T x 

then z T x ^ and 

a ■ z T Ayz T x 
a ■ z T Ay 

Thus p {AD(y)AD(z)) = a ■ \z T Ay\. Hence 

wBx{po(AD(y)AD(z)):y,ze{-l,i} n } 

= a ■ max ^\z T Ay\ : y,z G {— 1, l} n | = a ■ r(A). 

□ 

Now everything is set for Poljak and Rohn's reduction [5]. 

Theorem 4. Let A G M nxn be a non-singular matrix, let K be a positive real number and 
let A = (1/K) • J. Then r(A) > K if and only if the matrix interval [A -1 — A, A -1 + A] 
is singular. 

Proof. By Lemma [21 the matrix interval [A~ l — A,A _1 + A] is singular if and only if 
p (AD(y)AD(z)) > 1 for some y,z G {-1, l} n . By Lemma El p (AD(y)AD(z)) > 1 for 
some y,z G {—1, l} n if and only if r(A) > K . □ 

Corollary 5. RK1-MA TRIX-INTERVA L SINGULARITY is NP-hard. □ 

Remark. Poljak and Rohn [5] show that RK1-MATRIX-INTERVAL SINGULARITY 
belongs to the class NP by proving the existence of a singular matrix in every singular 
matrix interval, with a polynomial bound on the size of all entries of that matrix. 

4 Reduction from RK1-MATRIX-INTERVAL SINGULARITY 
to P-MATRIX 

The described reduction is by Coxson pQ. 



= Ax / 0, 



= A • z T x, 
= A. 
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Let A,Ag M. nxn . Consider the matrix interval [A, A + A]. Let A*'- 7 be the matrix 
whose element in the ith row and jth column is Ajj and which has zeros elsewhere. 
Then each matrix M in the interval [A, A + A] can be uniquely expressed as 

n 

M = A+ X>ijAH (1) 

i,i=l 

where pij G [0, 1] for all values of 

Each matrix A 1 '- 7 is a rank-1 matrix (even if A has higher rank), and so = fijsjj 
for some vectors rj Sj j € M n . We can actually take rj to be Aj in its ith entry and 
zero elsewhere, and Sjj to be 1 in its jth entry and zero elsewhere. 

Now let R be the matrix whose columns are all the n 2 vectors r^- and let S be the 
matrix whose columns are all the n 2 vectors Sij. Thus A = RS T . Moreover, if p G M" 2 
is the vector formed by the numbers pij, we can write ([I]) as 

M = A + RD{p)S J . 

Suppose that A is non-singular. Then the matrix interval [A, A + A] is non-singular 
if and only if 

det{A + RD(p)S T ) = det(A) det(J„ + A~~ 1 RD(p)S T ) ^ (2) 

for each vector p G [0, l] n . 

Supposing that the matrix A is non-singular, inequality ([2j) holds if and only if 

det(J„ + A -1 .RD(p)£ T ) / 0. (3) 

In this way we have proved that for a non-singular matrix A, singularity of the matrix 
interval [A, A + A] is equivalent to the existence of a vector p € [0, l] n that does not 
satisfy inequality Q. Since the expression in ([3]) is a multi-affine function of p, we can 
actually derive another condition. 

Lemma 6. Let ip(p) = det(/ n + A~ 1 RD(p)S T ). Then inequality (02) holds for each 

2 2 

p € [0, 1]™ if and only if ip(p) > for each p € {0, l} n . 

Proof. First observe that ip{p) = det(/„ + A~ 1 RD(p)S T ) is a multi-affine function of p, 
that is, for each i we have ip(p) = c\ + C2Pi, where ci, C2 depend on i and pj for j ^ i. 

We claim that any multi-affine function : [0, l] fc — > M. is non-zero on the whole 
domain if and only if its values on the vertices {0, l} k have all the same sign. Assuming 

2 

this claim holds, we notice that ip(0) = detl n = 1 > 0, so tp is non-zero on [0, 1]™ if and 
only if it is positive on {0, l} k . 

To prove the claim, first suppose that 4> is non-zero on [0, if but there are two vertices 
u, v G {0, l} k such that 4>{u) < and (j)[v) > 0. Following the path along the edges 
of {0,1}, we will find two vertices u',v' € {0,1} that differ in exactly one coordinate 
and such that (j>(u f ) < and <j)(v') > 0. Without loss of generality we may assume 
that u[ = and v[ = 1, while u[ = v\ for i > 2. Let x G [0, l] k be defined by 
x\ = <fi(u')/(cj>(u') — cj>(v')) and x\ = u\ for i > 2. Then cj){x) = 0, a contradiction. 

Conversely, if <fi is positive (negative) on all the vertices, it is easy to prove by induction 
on face dimension that <f> is positive (negative) in every internal point of each face. □ 
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Lemma [6] together with the discussion that precedes it imply the following character- 
isation. 

Lemma 7. Let A be a non-singular matrix and let R, S be defined as above. Then the 
matrix interval [A, A + A] is singular if and only if 

det(J n + A~ 1 RD(p)S T ) < 

for some p G {0, l} n . □ 

In order to get D(p) from the middle of the product to the beginning, we use the 
following lemma, whose proof we present in the Appendix. 

Lemma 8. Let F G R kxn and G G R nxk . Then det(I k + FG) = det(J„ + GF). □ 

This fact can be exploited to prove the following equivalence. 

Theorem 9. Let A be a non-singular matrix and let R, S be defined as in Lemma [7] 
Then the matrix interval [A, A + A] is singular if and only if the matrix M = I n 2 + 
S T A~ 1 R is not a P-matrix. 

Proof. Because of Lemma [8j 

ip(p) = det(/ n2 + A- 1 RD(p)S T ) = det(/ n 2 + D(p)S T A- 1 R). 

If p G {0, l}™ 2 and p ^ 0, the expression det(J n 2 + D{p)S J A~ l R) is equal to the 
principal minor of the matrix M obtained by selecting exactly those rows and columns 
that correspond to the 1-entries of the vector p. Thus ip(p) is non-positive for some 
p G {0, l} n if and only if the matrix M is not a P-matrix. 

The proof is now completed by applying Lemma □ 

Corollary 10. The problem P-MATRIX is co-NP- complete. 

Proof. NP-hardness follows from Corollary [5] and Theorem [9l 

The problem belongs to co-NP because after guessing the rows and columns, the 
corresponding principal minor, which certifies the negative answer, can be computed in 
polynomial time. □ 

Appendix: Proof of Lemma M 

One of the basic facts about determinants is that adding a multiple of a row to another 
row does not change the determinant. The following lemma (Theorem 3 in Section 2.5 of 
Gantmacher's book [3]) is a block version of this fact. Even though it holds for matrices 
with an arbitrary number of blocks, we state it just for 2x2 blocks. This variant is 
sufficient for the proof of Lemma [HJ 
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Lemma 11. Let A G W nxn be a matrix with block structure 



ni n 2 
mi{ / Ai : i A 1)2 
"*2 { A 2j2 



and let X G R m ^ xm ^ Y G M raiX?l2 . T/ien 



detA = detf^+^ A * + XA *A =detft' 1 

V A 2 ,i A 2j2 J \A 2 ,1 



Ai,i Ai )2 + Ai,iy 

a 2i2 + A 2 iy 



Proof. Since 

/Ai,! + XA 2>1 A lj2 + XA 2j2 \ // mi I\ 

v am a 2j2 y v o w 

we have 



det( A lil +XA 2>1 A li2 + XA 2 , 2 \ =de // X\ 

A 2 ,i a 2i2 y \ 



Similarly 



dr( A^ + A^yX = // ^j _ (lr( | n 

\A 2il a 2i2 + a 2i1 y-> I / ' 



Finally comes the proof of Lemma 
Proof of LemmalM Applying Lemma [TT] twice, we get 

det(4 + FG) = det( /fc + FG ?)-^(n /) 

" d6t (g J n + Gf) = det(/n + GF) - 

Here (*) follows by applying Lemma PTT1 to rows with X = F and (f) follows by applying 
it to columns with Y = F. □ 
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